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ABSTRACT 

Context. M31N 2015-01a (or M31LRN 2015) is a red nova that erupted in January 2015 - the first event of this kind observed in M31 
since 1988. Very few similar events have been confirmed as of 2015. Most of them are considered to be products of stellar mergers. 
Aims. Results of an extensive optical monitoring of the transient in the period January-March 2015 are presented. 

Methods. Eight optical telescopes were used for imaging. Spectra were obtained on BTA, GTC and the Rozhen 2m telescope. 
Results. We present a highly accurate 70 d lightcurve and astrometry with a 0.05"uncertainty. The color indices reached a minimum 
2-3 d before peak brightness and rapidly increased afterwards. The spectral type changed from F5I to F0I in 6 d before the maximum 
and then to K3I in the next 30 d. The luminosity of the transient was estimated to 8.7 4 2'] x 10 5 L o during the optical maximum. 
Conclusions. Both the photometric and the spectroscopic results confirm that the object is a red nova, similar to V838 Monocerotis. 

Key words, stars: novae, cataclysmic variables - stars: individual (M31N 2015-01a, M31-RV, V838 Mon, V1309 Sco, V4332 Sgr) 


1. Introduction 


Red novae (or red transients) are a rare class of optical tran¬ 
sients, reaching a peak luminosity equal to or higher than the 
brightest classical novae, but still lower than supernovae. Many 
of t h ern are best explained by binary star mergers ( Tv lend a & 
Soker 2006) and therefo re are often referred to as stellar merg - 
ers or "mergebursts" (iTvlenda et al.ll201 ll : lKochanek et al .120141) . 
Common observational properties during the first months of the 
outburst are the initial slow decline and a spectrum changing to¬ 
wards a late-type supergiant phase with colors shifting to the red. 

Modern observational history of these events goes back to 
1988 when the transient M31-RV was detected bv [Rich et alJ 
( 1989 ) as an M0 supergiant variable with peak absolute M\ w \ ~ 
-10. A similar outburst in the Milky Way was observed in 1994 
- V4332 Sgr changed its s pectrum from K3 to M8-9 during a 3- 
month period ( Martini et al.l 19991) . The eruption of V838 Mon in 
February 2002 (iMunari et al.l2002l) reached an apparent magni¬ 
tude of V=6.7 mag and was studied quite closely. The lightcurve 


showed three maxima and a consistent change of colors red- 
wards as the spectrum reached an L supergiant phase at the end 
of the year (lEvans et al.ll2003i) . The resulting circumstellar light 
echo was used to obtain a distance estimate of 6.1±0.6kpc, cor¬ 
responding to a peak abs olute magnitude My = -9.8 ( Bon d et al.l 

2003: [Sparks et al 120081) . A bright (peak M r -12) event, M85 

OT2006-1, was reported bv lKufkarni et al.l (120071) . Like other red 
novae, it showed a narrow Her emission, a temperature decrease 
and an infrared excess a few months after the eruption (iRau et ali 
120071) . 

VI309 Sco is another Milky Way red nova, discovered in 
20Q8. lTVlenda et all ( 2011 ) have shown that the progenitor is a 
contact binary with a 1.4-day period, thereby adding a significant 
argument in favour of the stellar merger model for the cause of 
these outbursts. 

It has been suggested thatnhe first catalogued nova, CK Vul 
was also a "red transient" (Kaminski et al.ll2015 ). A p ossible red 
nova in M101 was reported bv lGoranskii et"alJ ( 201 5!). 
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Table 1. Telescopes used for the monitoring of M3 IN 2015-01a. The 
last column contains the designations used in the online Table[4] 


Table 2. List of reference stars. The magnitudes are calibrated by the 
iMassev et al. (12006 ) catalog using images from the 6m BTA, the Rozhen 
2m and the Konkoly Schmidt telescopes. 


Telescope 

Observatory 

Country 

Image scale 

Des. 

a (J2000) 

8 (J2000) 

B 

V 

R 

I 




[arcsec/px] 


00:42:08.9 

+40:55:33 

14.521 

13.739 

13.312 

12.886 

6 m BTA 

SAO RAS 

Russia 

0.36 

BTA 

00:42:12.6 

+40:54:38 

14.792 

14.267 

13.996 

13.699 

2.5 m NOT 

ORM 

Spain 

0.19 

NOT 

00:42:12.2 

+40:55:49 

16.030 

15.287 

14.883 

14.495 

2 m RCC 

Rozhen NAO 

Bulgaria 

0.74 

R2M 

00:42:05.3 

+40:53:38 

16.876 

15.915 

15.381 

14.822 

1 mLCOGT 

McDonald 

USA 

0.47 

LCOGT 

rnrs of calibration: 


0.015 

0.007 

0.008 

0.015 

65 cm 

Ondrejov 

Czech Rep. 

1.05 

ONDR 







60/90 cm Schmidt 

Konkoly 

Hungary 

1.03 

KSCH 







50/70 cm Schmidt 

Rozhen NAO 

Bulgaria 

1.08 

RSCH 







50 cm Cassegrain 

Uni Leicester 

UK 

0.89 

UL50 

-'-'- 

—I ->- 

— - - r 

-t- 

- T - 



Here, we present a new addition to this class of objects, a red 
nova that erupted in M31 in Jan 2015. M31N 2015-01a was dis¬ 
covered by the MASTER-Kislovodsk auto-detection system on 
2015 Jan 13 at ~19.0mag, unfiltered (Shurnkov_et_aL 2015). It 
was in correctly identified as a classical nova bv lKurtenkov et alJ 
(12015al) based on a strong Ha emission on the Jan 16 spectrum, 
also reported by iHodgkin et al.1 (i2015l) . On Jan 15 the object 
sho wed the spectral charact eristics of an F5 s upergiant ( Fabrika 


et al. l2015h l The transient became very bright, reaching a peak 
magnitude of R~ 15.1 mag circa Jan 22. It then proceeded to fade 
slowly and consistently turned redder. At the end of February the 
colors had changed considerably with the spectrum resembling 
that of a K supergiant. Subsequently, we announced that the tran¬ 
sient is a red nova (iKurtenkov et al.ll2015bt) . Our supposition that 
the overall spectral evolution of the object is very similar to the 
observed in V838 Mon was confirmed by [Williams et al.1 d2015l 
Fig. 2). The poss ible progenitor system ha s been discussed by 
iDong et aD (1201 5h and lWilliams et al.1 (120151) . 
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Fig. 1. Lightcurve of M3 IN 2015-01a in Johnson-Cousins BVRI filters. 
The zero-point moment roughly coincides with the R-band peak. The 
absolute magnitude scale (right) applies to the R-band datapoints only. 


2. Observations and data reduction 

We performed imaging of M3 IN 2015-01 a in Johnson-Cousins 
BVRI filters (Table|T|). The science frames were divided by flat 
fields after dark/bias subtraction. 

Images of vastly different depths and FOVs were obtained. In 
order to do all the photometry in a similar manner, we used four 
bright reference stars, within 1.6' from the transient (Tabled. 
That way systematic effects over large FOVs are not affecting 
the results. This also allowed us to minimize the instrumental er¬ 
rors of the references. Three out of those four stars were used for 
each frame. [Massev et al . d2006 t ) found systematic differences of 
the order of 0. 1 mag between their magnitudes and the ones by 
iMagnier et al.1 (119921) . so we used the Massey catalog to recali¬ 
brate the BVRI magnitudes of the four references. Two images 
per filter and between 18 and 30 stars per image were used for 
the recalibration. The peak signals of the four reference stars on 
those images are within the linear range of the detector. The re¬ 
sults are presented in Table[2] 

Aperture photometry was done for all individual frames. The 
aperture radii we used were in the range of 1-1.5 FWHM of the 
Gaussian profiles of point sources. The results are presented as 
online material in Table[4] Only the mean JD and median mag¬ 
nitude are given where several consecutive images are available. 

Spectra of M3 IN 2015-01 a were obtained with three tele¬ 
scopes - the 6 m BTA at SAO RAS, the 10.4 m GTC at ORM 
and the 2 m RCC telescope at Rozhen NAO (Table[3]). All the 
data reduction and calibrations were carried out with standard 


Tabl e 3. Log of spectral observations. The instrum ents are described by 
iJockers et all 1200(f) . (Afanasiev & Moiseev 120051) and lCepal 1 1 9981) . 


Date 

Telescope/Instrument 

R 

Range [A] 

Exp. time [s] 

Jan 15 

BTA/SCORPIO 

1000 

4050-5850 

1500 

Jan 16 

2mRCC/FoReRo2 

400 

5500-7500 

4x900 

Jan 17 

2mRCC/FoReRo2 

400 

5500-7500 

5x900 

Jan 21 

BTA/SCORPIO 

1000 

4050-5850 

1800 

Feb 21 

BTA/SCORPIO 

1000 

3720-5530 

2700 

Feb 22 

BTA/SCORPIO 

1000 

5750-7500 

1800 

Feb 24 

GTC/OSIRIS 

2500 

4430-9090 

6x180 


MiD.vjj] (BTA) and irai0 (2m RCC and GTC) procedures. The 
BTA and the GTC spectra were calibrated in relative fluxes using 
spectrophotometric standards G191-B 2B, Hiltner 600 and Hz 2 
(lOkdl 1 9901: lHamuv et akl[l99l[l994l) . 

3. Results and discussion 

3.1. Astrometry 

Astrometric calibra tions were made using both the PP MXL 
dRoeser et alJl2010h and the LGGS (IMassev et all [2006) cata- 

1 See http://www.eso.org/sci/software/esomidas/ 

2 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy (AURA) under a cooperative agreement with the National 
Science Foundation. 
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Fig. 2. The change of the observed color indices. Vertical lines represent 
the times the BTA (thin) and GTC (thick) spectra were taken at. 


logs. Astrometric solutions were obtained for 5 images from the 
10.4 m GTC, 6 m BTA and 1 m LCOGT telescopes with an 
atmospheric seeing in the range of 1.3" - 1.7". Separate cal¬ 
ibrations were made using 20 relatively bright stars from the 
PPMXL catalog and 16 fainter stars close to the transient from 
LGGS. From PPMXL we derived a, <5(J2000) = 0() h 42 m 08:053, 
+40°55'01727 with uncertainties of Aq. ~ 0.2", A s ~ 0.1". 
The LGGS calibrations yielded a,<i(J2000) = 00 h 42 m 08?065, 
+40°55'01'.'33 with A* ~ A* ~ 0.05". 

3.2. Photometry 

The light curve of M3 IN 2015-01 a (Fig. Q} shows a distinct 
similarity to the VI309 Sco lightcurve during the first month of 
the outburst dMason et al.ll2010h . The Jan 21 datapoints are ob¬ 
tained shortly after the B-band and before the /Gbaud maxima, 
so we assume that they coincide with the V-band maximum at 
15.43 mag. The initial fading has been quite slow. The decline 
times by 2 mag after the maximum ti were ~ 17 c/, ~ 40 d and 
~ 50 d in B. V and /Gbauds respectively. The colors have been 
consistently shifting to the red, e.g. B - R increased by 2.1 mag 
in the 31 days after the maximum. The R magnitude reached a 
plateau 0.6mag below the maximum in ~ 10d. but no rebright¬ 
ening followed, as was the case with V838 Mon. 

Figure[2] shows the change of colors. The color temperature 
initially increased and reached a peak 2-3 d before the /Gband 
maximum. Afterwards it has consistently decreased as the B - R 
index increased by ~2 mag in one month. 

3.3. Spectroscopy 

The aim of the FoReRo2 spectra of M3 IN 2015-01 a obtained 5 
and 4 days before the maximum was to confirm that the object 
is a nova. The only features in these spectra are the emission 
lines of Her with equivalent widths 19 + 1 A and 17 + 1 A on 
Jan 16 and 17 respectively and a Nai doublet absorption. After 
a car eful re-calibration of the wavelength scales (see Kurtenkov 
et al. 12015ah the measured radial velocities of the Her emission 
peaks are -365 km s 1 and -350km s' 1 . In good agreement with 
the velocity of ~ -370 ± 50 km s' 1 estimated by iFabrika et al.l 
( 2015 1 for the part of M31 where the object is located. 

The first SCORPIO spectrum was obtained 6 days before 
the maximum brightness. The spectral region covered (4050- 
5850 A) is dominated by absorption lines of ionized and neu¬ 
tral elements like Fe n, Ti n, Cru, Mg n and Fe i. Hy and He) are 


visible in absorption, while H j8 absorption is partially filled by 
an emission component. A comparison of the absorption line 
spectrum with the iJacobv et al.l d 19841) library of stellar spec¬ 
tra shows very close similarity with a F5I spectrum. The next 
SCORPIO spectrum obtained on Jan 21 more or less coincides 
with the brightness maximum. The absorption line spectrum re¬ 
mained practically the same. However, the continuum is much 
bluer (Fig. [3]) showing a better resemblance with a F0I spec¬ 
trum. We cross-correlated these two spectra with a F5I tem¬ 
plate and found a radial velocity for the pre-maximum spectrum 
-462 + 38 km s' 1 and -534+14 km s' 1 for the spectrum around 
brightness maximum (measured using the iraf fxcor package). 
The cross-correlation of the M3 IN 2015-01a spectra itself, con¬ 
firms a velocity difference of 64 ± 13 km s 1 . Taking into consid¬ 
eration the velocity -370 km s' 1 , connected with M31, we found 
that the heliocentric velocity of the ejected matter on Jan 15 was 
~ 90 km s' 1 and increased to ~ 160 km s' 1 on Jan 21. 

One month after the maximum we obtained two additional 
SCORPIO spectra in consecutive nights and one OSIRIS spec¬ 
trum two days later. At that time the V brightness of the star 
decreased by ~ 1.5 mag and the colors drastically reddened. Cor¬ 
respondingly, the spectrum of M3 IN 2015-01 a moved to the 
lower temperature classes. Now, the neutral elements absorp¬ 
tions dominate the spectrum together with weak TiO molecular 
bands. Among the strongest absorptions in the spectrum are the 
Nai doublet. Ban 6142A and 6497A Cai 6573 A confirmed 
bv lWilliams et ali (120151) . and the near IR triplet of Can. The 
Can H & K lines are detected in emission. The strong Nai lines 
are well split in our OSIRIS spectrum. This allowed us to esti¬ 
mate their FWHM~ 200 ±15 km s' and their radial velocity to 
-570 ± 15 km s' 1 . In good accordance with the velocity mea¬ 
sured in all our February spectra of M31N 2015-01a (see be¬ 
low) and testifies to the stellar origin of these absorptions. In the 
OSIRIS spectrum Her presents as a very weak emission divided 
in two components by a central absorption. In the SCORPIO 
spectrum only a marginal emission component can be seen. A 
comparison of the M31N 2015-01a spectrum about one month 
af ter the brightn ess maximu m with the stel l ar sp ectra libraries 
of IJacobv etalJ (1 19841) and lLe Borgne et all (120031) shows clos¬ 
est similarity with a K3-4I spectral type. Cross-correlating the 
February spectra with a K3I template we estimated an aver¬ 
age velocity -588 + 40 km s' 1 , suggesting an increase of the 
ejected matter velocity to ~220kms '. The cross-correlation of 
the SCORPIO spectra with the OSIRIS one shows small velocity 
differences of 8.4 ± 6.6 km s 1 and -10 ± 14 km s 1 for Feb 21 
and 22 respectively. All uncertainties are ler equivalent. 

Using the task fitspec in the stsdas_ synphot0 and tem- 
plate s pect ra from lJacoby e t alJ ( 1984), Sanchez-Blazauez et al.l 
( 2006 ) and lLe Borgne et ak 120031) we fitted our spectra in an at¬ 
tempt to estimate the reddening. Fitting the spectrum on Jan 15 
with a F5I template we obtained a reddening £b-v = 0.37 mag. 
For the spectrum on Jan 21 fitted with a F0I template the corre¬ 
sponding £b-v was 0.33 mag. Merged spectra on Feb 21 and 22 
and the spectrum on Feb 24 were best fitted with a K3I template 
and Eb-v = 0.25 mag and Eb-v = 0.42 mag respectively. 

Combining the column hydrogen densities of Nieten et al.l 


(2006) with the gas-to-dust (fNedialko v et alJl201 1 ) ratio yields 
total color excess on the line of sight Eb-v = 0.44 + 0.08 mag. 
An identical result: Eb-v = 0.42 ± 0.03 mag is derived from the 
maps of the M31 dust surface density riDraine et alil2014l) . A 
comparison of the B — V and V — R colors with the catalog com- 


3 stsdas is a product of the Space Telescope Science Institute, which 
is operated by AURA for NASA. 
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Fig. 3. BTA/SCORPIO spectra obtained with the VPHG1200B and G 
grisms. The Jan 21 spectrum is close to maximum brightness and arbi¬ 
trary shifted to unity. The relative flux and shape of the Jan 15 and Feb 
21 spectra reflect the evolution of the luminosity and the colors. 
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Fig. 4. The combined GTC/OSIRIS spectrum dereddened by E B -v = 
0.4 mag in comparison to HD 157999 (K3Iab) and 41 Gem (K4Iab). 


piled by [Worthev & Lee! (1201 ll) shows that the reddening Eb-v 
is in the range of 0.15 - 0.58 mag. The mean reddening value, 
obtained from the spectra - Eb-v = 0.35 + O.lOmag, of which 
-0.05 mag can be attributed to the interstellar reddening in our 
Galaxy dSchlaflv & Finkbeinerl20l7t) . is well within this range. 
It suggests a true color of ( B - V)o = 0.15 + O.lOmag and a V- 
band bolometric correction of +0.03 mag at maximum. Adopt¬ 
ing a distance modulus (m - M) = 24.47 ma^Q and a total-to- 
selective extinction Ry = 3.1 yields a peak absolute magnitude 
My = -10.13 ±0.30mag. This sets the luminosity limits at max¬ 
imum on Lboi = 8.7+2 2 x 1()5 l o- Note, that M3 I N 2015-01a is 
projected inside a tiny (2"x2") HII region #1527 (lAzimlu et al,l 
1201 lh . The latter list an extinction A R = 0.538 mag, which corre¬ 
sponds to the somewhat lower Eb-v = 0.23 mag. An association 
with this HII region could decrease the luminosity estimate to 
~ 6 x 10 4 5 L o . 


4. Summary 

The outburst of M3 IN 2015-01 a was a very luminous event, 
reaching ~ 10 6 L o at the optical maximum. We present accurate 
photometry and astrometry in the period January-March 2015. 

4 The NED database (http: //ned. ipac. caltech. edu/1 lists a me¬ 

dian (m - M) = 24.45 mag for M31. If the object lies in the M31 plane, 

its location suggests a +0.02 mag difference from the center. 


The transient shares all major observational characteristics of red 
novae and we conclude that it belongs to this group of events. 
The fading was slow at first with D - 50 d in /Gbaud. In a one- 
month period B -R increased by 2.1 mag and the observed spec¬ 
tral type changed from F0I to K3I. The exploration of red novae 
is so far based on a very small sample of sources, which makes 
M31N 2015-01a so important: it moves us a step closer to rec¬ 
ognizing similarities and differences in their behaviour. 
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Table 4. Calculated magnitudes for M31N 2015-01a in Johnson- 
Cousins BVRI filters. 


JD-2457000 filter 

mag 

err 

telescope 

airmass 

36.264 

B 

19.641 

0.161 

RSCH 

1.17 

38.248 

B 

17.046 

0.017 

BTA 

1.30 

38.345 

B 

16.992 

0.029 

RSCH 

1.69 

39.261 

B 

16.549 

0.017 

RSCH 

1.19 

39.413 

B 

16.498 

0.029 

UL50 

1.53 

40.295 

B 

16.215 

0.017 

RSCH 

1.35 

44.196 

B 

15.934 

0.017 

BTA 

1.16 

46.612 

B 

16.109 

0.026 

LCOGT 

1.32 

48.658 

B 

16.679 

0.242 

LCOGT 

1.71 

51.611 

B 

17.261 

0.026 

LCOGT 

1.38 

60.323 

B 

17.865 

0.096 

KSCH 

1.79 

61.252 

B 

17.974 

0.026 

KSCH 

1.32 

67.257 

B 

18.294 

0.066 

KSCH 

1.42 

68.246 

B 

18.383 

0.022 

RSCH 

1.53 

72.249 

B 

18.719 

0.103 

KSCH 

1.46 

75.224 

B 

18.791 

0.018 

BTA 

1.92 

93.340 

B 

20.076 

0.051 

NOT 

2.53 

96.346 

B 

20.385 

0.041 

NOT 

2.99 

38.246 

V 

16.416 

0.008 

BTA 

1.29 

44.187 

V 

15.436 

0.011 

BTA 

1.14 

66.319 

V 

16.626 

0.012 

KSCH 

1.93 

67.243 

V 

16.621 

0.026 

KSCH 

1.35 

72.241 

V 

16.832 

0.065 

KSCH 

1.41 

75.225 

V 

16.873 

0.013 

BTA 

1.93 

89.323 

V 

17.869 

0.059 

UL50 

1.92 

92.246 

V 

18.098 

0.108 

KSCH 

1.91 

36.243 

R 

18.087 

0.026 

RSCH 

1.11 

37.210 

R 

17.213 

0.011 

R2M 

1.05 

38.246 

R 

16.051 

0.011 

BTA 

1.29 

38.330 

R 

16.047 

0.009 

RSCH 

1.55 

39.242 

R 

15.679 

0.009 

RSCH 

1.13 

40.283 

R 

15.407 

0.009 

RSCH 

1.29 

40.390 

R 

15.364 

0.012 

R2M 

2.38 

44.189 

R 

15.133 

0.011 

BTA 

1.14 

46.617 

R 

15.165 

0.013 

LCOGT 

1.34 

48.661 

R 

15.347 

0.032 

LCOGT 

1.75 

51.615 

R 

15.609 

0.011 

LCOGT 

1.41 

57.285 

R 

15.678 

0.011 

ONDR 

1.36 

59.238 

R 

15.692 

0.011 

ONDR 

1.20 

60.317 

R 

15.675 

0.012 

KSCH 

1.73 

61.234 

R 

15.683 

0.011 

KSCH 

1.24 

62.596 

R 

15.681 

0.009 

LCOGT 

1.50 

64.608 

R 

15.689 

0.009 

LCOGT 

1.66 

66.310 

R 

15.725 

0.011 

KSCH 

1.82 

67.239 

R 

15.741 

0.009 

KSCH 

1.33 

68.230 

R 

15.804 

0.013 

RSCH 

1.42 

72.244 

R 

15.839 

0.027 

KSCH 

1.42 

73.250 

R 

15.858 

0.009 

ONDR 

1.40 

75.223 

R 

15.894 

0.024 

BTA 

1.90 

80.270 

R 

16.111 

0.011 

ONDR 

1.68 

86.329 

R 

16.417 

0.021 

ONDR 

2.72 

87.302 

R 

16.437 

0.016 

ONDR 

2.28 

89.266 

R 

16.642 

0.012 

ONDR 

1.88 

89.337 

R 

16.617 

0.061 

UL50 

2.08 

92.245 

R 

16.987 

0.104 

KSCH 

1.90 

93.345 

R 

17.011 

0.011 

NOT 

2.68 

96.351 

R 

17.315 

0.013 

NOT 

3.19 

99.283 

R 

17.617 

0.035 

ONDR 

2.54 

105.279 

R 

18.499 

0.084 

ONDR 

2.79 

60.320 

1 

15.028 

0.021 

KSCH 

1.76 

61.245 

I 

15.017 

0.018 

KSCH 

1.28 

67.241 

I 

15.019 

0.018 

KSCH 

1.34 

68.272 

I 

15.085 

0.016 

RSCH 

1.78 

72.246 

I 

15.048 

0.028 

KSCH 

1.44 

89.357 

I 

15.344 

0.111 

UL50 

2.34 

92.248 

I 

15.377 

0.065 

KSCH 

1.92 

93.350 

I 

15.557 

0.021 

NOT 

2.84 

96.356 

I 

15.735 

0.021 

NOT 

3.42 



